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Online monitoring of viable cell volume (VCV) is essential to the development, monitoring, 
and control of bioprocesses. The commercial availability of steam- sterilizable dielectric- 
spectroscopy probes has enabled successful adoption of this technology as a key noninvasive 
method to measure VCV for cell-culture processes. Technological challenges still exist, how- 
ever. For some cell lines, the technique s accuracy in predicting the VCV from probe- 
permittivity measurements declines as the viability of the cell culture decreases. To investi- 
gate the cause of this decrease in accuracy, divergences in predicted vs. actual VCV meas- 
urements were directly related to the shape of dielectric frequency scans collected during a 
cell culture. The changes in the shape of the beta dispersion, which are associated with 
changes in cell state, are quantified by applying a novel "area ratio" (AR) metric to 
frequency-scanning data from the dielectric-spectroscopy probes. The AR metric is then used 
to relate the shape of the beta dispersion to single-frequency permittivity measurements to 
accurately predict the offline VCV throughout an entire fed-batch run, regardless of cell 
state. This work demonstrates the possible feasibility of quantifying the shape of the beta dis- 
persion, determined from frequency-scanning data, for enhanced measurement of VCV in 
mammalian cell cultures by applying a novel shape-characterization technique. In addition, 
this work demonstrates the utility of using changes in the shape of the beta dispersion to 
quantify cell health. © 2013 American Institute of Chemical Engineers Biotechnol. Prog., 
30:479-487, 2014 
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Introduction 

Online monitoring is essential to the development, monitor- 
ing, and control of bioprocesses. 1 ^ In bioreactors, one of the 
most important online measurements is total biomass, or via- 
ble cell volume (VCV), because concentration measurements 
of the media must be normalized to VCV to yield useful 
information about the cell state. 5-9 Several methods have been 
developed in an attempt to achieve this observability online, 
including oxygen-demand analysis, 10 turbidity analysis, 11 and 
dielectric spectroscopy, 12-14 as well as at-line methods (i.e., 
sterile autosampling) and more-traditional methods (e.g., try- 
pan blue dye exclusion). 15 Of these methods, dielectric spec- 
troscopy offers numerous advantages including sensitivity to 
cell numbers, a high scan rate, and maintenance of bioreactor 
sterility, as well as the capability to measure cell properties 
not easily accessed by other online methods. 12,16 

In dielectric spectroscopy, the dielectric permittivity of 
cells in a population — that is, their ability to store electrical 



charge — is measured as a function of the frequency of the 
applied electrical field using a dielectric probe. Any cell 
with an intact membrane can act as a dielectric material 
when subjected to an electric field by means of a well- 
known phenomenon called Maxwell-Wagner polarization. 17 
Since each cell stores electrical charge independently at cell 
densities relevant to a bioreactor, the total permittivity is 
proportional to the cell volume fraction. 

In industrial settings, dielectric probes are mainly used to 
measure total VCV online in bioreactors. For this applica- 
tion, permittivity is usually measured at only one or two fre- 
quencies. However, dielectric spectroscopy has also been 
investigated extensively in academic settings to measure the 
permittivity characteristics of dielectrics across a frequency 
spectrum using frequency scanning. Several informative 
reviews have been published on the behavior of the dielec- 
tric spectrum as it pertains to cells. 17-20 Some of the key 
physics are reviewed briefly here. 



This is an open access article under the terms of the Creative 
Commons Attribution-NonCommercial-NoDerivs License, which permits 
use and distribution in any medium, provided the original work is 
properly cited, the use is non-commercial and no modifications or 
adaptations are made. 

Correspondence concerning this article should be addressed to 
Brandon J. Downey at brandon.downey@bendresearch.com. 



Background 

Overview of frequency scanning 

If the frequency of the applied electric field is varied, a 
frequency-dependent response can be observed in the 
measured permittivity values. This characteristic response 
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Figure 1. Representation of the beta dispersion. 



(or dielectric relaxation) is known as the beta dispersion, and 
occurs mainly between 0.1 and 100 MHz for mammalian cells. 
Beta dispersions may consist of a single relaxation or a sum 
of multiple relaxations, as illustrated in Figure 1. Dielectric 
relaxation phenomena are commonly characterized by a Cole- 
Cole model, in which dielectric relaxations are characterized 
by a dielectric increment (Ae), critical frequency (f c ), and a 
Cole-Cole parameter (a). 21 

The dielectric increment is defined as the difference 
between the low- and high-frequency plateaus of the meas- 
ured beta dispersion. The critical frequency is defined as the 
frequency at which the inflection occurs in the beta disper- 
sion and is directly related to the average dielectric relaxa- 
tion time of the sample (i.e., f c = 2m, where x is the 
dielectric relaxation time). The Cole-Cole parameter is an 
empirical parameter that determines the "steepness" with 
which the permittivity decreases with frequency and is 
thought to be related to the distribution of dielectric proper- 



ties of the measured sample, 
introduced in Eq. 1: 



The Cole-Cole equation is 



A £ [l + (/// c ) (1 " a) cos(7r(l-a)/2) 
l + 2(^// c )( 1 - a )cos(7r(l-a)/2) + (f// c ) 2 ( 1 - a ) ! 



(1) 



where s(f) is the permittivity as a function of frequency if) 
and e h is a permittivity offset equal to the value of the high- 
frequency plateau of the beta dispersion. The values of the 
parameters of the Cole-Cole model describing the shape of 
the beta dispersion depend on the properties of the dielectric 
being measured. Previously, many cell properties have been 
observed to affect the shape of the beta dispersion, including 
cell size 23 ; intracellular organelle content 24 ' 25 ; membrane- 
specific capacitance ' ; and intracellular conductivity. ' 
Changes in the shape of the beta dispersion have also been 
reported for healthy versus nonhealthy cells 28-32 and for cells 
undergoing other metabolic changes. ' ' 



Overview of VCV measurement techniques using dielectric 
spectroscopy 

Information about the shape of the beta dispersion can be 
obtained by measuring capacitance as a function of fre- 
quency at many points along the beta dispersion, instead of 
measuring at only one or two frequencies. However, for sim- 
plicity, VCV is often estimated by measuring the beta dis- 
persion only at one or two frequencies. 34 When only two 
frequencies are used, the frequencies can be chosen so that 



the resulting permittivity values at those frequencies are 
indicative of the dielectric increment. 

The dielectric increment can also be estimated by meas- 
uring at a single frequency if the high-frequency plateau of 
the beta dispersion is assumed to be constant. The dielectric 
increment is known to relate directly to the total bio volume, 
or cell volume fraction, of the system via the relationship. 34 



Ae = 



(2) 



where r is the average cell radius, <P is the total cell volume 
fraction, and C m is the membrane-specific capacitance (usu- 
ally assumed to be a constant). 

VCV, therefore, is usually predicted assuming that (1) a 
linear correlation exists between VCV and dielectric incre- 
ment and (2) all remaining terms in the Cole-Cole model 
describing the beta dispersion are constant over the entire 



experiment, 34 as shown in Eq. 3: 

VCV (f)=AXAe(f)+*i, 



(3) 



where VCV(0 is the predicted viable cell volume at a given 
time (f), A is a constant of proportionality relating the VCV 
to Ae, and k x is a constant offset between VCV and Ae. 

The commercial availability of steam-sterilizable dielectric 
probes has allowed dielectric spectroscopy to be successfully 
used in a number of industrial applications for estimating 
VCV and for determining the timing of feeding and tempera- 
ture shifts in bioreactors. Despite these successes, the full 
value of dielectric spectroscopy has not yet been realized for 
use in applications such as continuous feed control. 16 One rea- 
son may be the large errors that have been observed for some 
cell lines when trying to predict VCV with this method. 23 ' 35 ' 36 

Specifically, for some mammalian cell cultures, the single- 
frequency permittivity correlates linearly during the growth 
phase with the results of offline viable cell-counting methods 
based on trypan blue dye exclusion. As the viability of the 
culture declines, however, the linear correlation (Eq. 3) 
becomes progressively worse at accurately predicting the 
actual VCV. This phenomenon has been observed previously 
with dielectric probes monitoring Chinese hamster ovary 
(CHO) cells. 23 ' 35 ' 36 In most published cases, the measured 
permittivity overpredicts the actual VCV in the death phase 
of the culture. In at least one case, however, only negligible 
amounts of divergence were observed in the death phase. 37 

By linearly correlating the dielectric increment to VCV, 
the dielectric properties of each cell is assumed to be con- 
stant over time. In cell-culture processes, particularly in fed- 
batch mode, cell properties are known to change, sometimes 
drastically, over the course of the culture as nutrients are 
depleted, metabolic by-products accumulate, and cell aging 
occurs. Changes in cell size, 38 metabolic state, 39 viability, 
and growth rate, among other cell properties, have all been 
observed to change over the course of fed-batch cell culture. 
Some of these cell-property changes — particularly cell 
size, ' metabolic state, and degree of apoptosis ' — 
have been linked to shifts in the dielectric properties of the 
cell population. If the dielectric properties of the cells 
change significantly over time or if nonviable cells produce 
a dielectric signal, then a divergence would be expected as 
the cells enter the death phase of the culture, due to the 
changing dielectric properties of the cells. The divergence 
phenomenon may be caused by nonviable cells contributing 
additional permittivity signal to the measured dielectric 
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increment (i.e., cells that are counted as nonviable by trypan 
blue dye exclusion are still participating in Maxwell-Wagner 
polarization). 

In the current work, it is hypothesized that nonviable cells 
do produce a permittivity signal. It is also hypothesized that 
the permittivity signal produced by nonviable cells is distin- 
guishable from the signal produced by viable cells in terms 
of the frequency dependence of the signal (i.e., the shape of 
dielectric spectrum). If this is true, a correction is required to 
compensate for this additional permittivity signal to obtain a 
permittivity measurement that accurately predicts VCV and 
matches the results of the offline method. This article 
describes a novel method that accomplishes this. 

Materials and Methods 

Cell culture 

A GS — /— NS/0 cell line expressing a recombinant pro- 
tein was grown in a fed-batch culture in 2-L glass, autoclav- 
able bioreactors (Broadley- James Inc., Irvine, CA). The 
cultivation was inoculated at 2.5 X 10 5 cells/mL. Operating 
conditions were as follows: temperature 36.5°C ± 0.1°C, pH 
7.2 ±0.1, dissolved oxygen 30% ± 5% air saturation, and 
agitation 170 ± 5 rpm. A one-pitched blade impeller with 
three 45° blades was used (Metrohm Applikon BV, Schie- 
dam, The Netherlands). The culture was fed with a chemi- 
cally defined nutrient feed at a continuous flow rate after 
Day 3. Glucose was continuously controlled separately at 
1 ± 0.5 g/L. 

Offline cell viability and VCV measurements 

Sampling was performed at least once daily. Samples were 
analyzed on a Nova Bioprofile 400 (Nova Biomedical, Wal- 
tham, MA; data not shown) and a Cedex automated cell coun- 
ter (Roche Diagnostics Co., Indianapolis, IN). Samples were 
retained for additional analysis of metabolite concentrations 
and product titer (data not shown). Cell viability was deter- 
mined with the Cedex instrument by trypan blue dye exclu- 
sion. Cedex images were visually spot-checked for counting 
and edge-finding accuracy, and determined to be sufficiently 
accurate. The VCV was determined from the Cedex viable 
cell density (VCD) and average viable cell diameter (D) data 
approximated as spheres using 

4 (D\ 3 

VCV = -7r(-J XVCD, (4) 

where n is the numerical constant. Calculating VCV using 
visual cell-counting methods inherently assumes a spherical 
cell morphology, which is not necessarily a valid assumption 
for mammalian cells at low viabilities. This method for deter- 
mining VCV is based on using an average viable cell diame- 
ter calculated from a complex distribution of measured cell 
diameters. It is known that cell size distributions change 
dynamically over the course of a cell culture process 38 and 
changes in cell size distribution can be observed in the data 
collected here as well. In general, the size distribution 
appears to deviate significantly from a normal distribution at 
all points during the culture, with the deviation becoming 
even more pronounced during the death phase. Although 
using the mean of a complex distribution of cell sizes to cal- 
culate VCV is not ideal, gross changes in the average size 
distribution are still captured using this method, allowing 



demonstration of this proof of concept of the relationship 
between VCV and dielectric property measurements. 

Other methods may yield more-accurate measurements of 
VCV. An example of another useful offline measure that is 
indicative of VCV is viable packed-cell volume (vPCV). 35 
vPCV does not require image analysis to explicitly quantify 
the cell morphology for determination of cell volumes, so it 
does not require cell morphology or a distribution of diame- 
ters to be assumed. However, it is affected by the dynamics 
of cell packing as well as the bio volume. 

Collection of flow cytometry data 

Flow cytometry data were collected on a fluorescence- 
activated cell sorter (FACS)-Aria instrument (BD Bioscien- 
ces, San Jose, CA). Forward-scatter and side-scatter data 
were collected and used to gate cells into populations of 
"healthy" and "unhealthy" cells. The gated subpopulations 
had distinct side-scatter properties when gated according to 
the method described by Geerts et al. 41 As previously 
observed, 41 the "healthy" population linearly correlated to 
viable, nonstaining cells observed with trypan blue dye exclu- 
sion (R 2 = 0.95). The "unhealthy" population correlated to 
cells that stain positive (i.e., nonviable) in trypan blue dye 
exclusion (R 2 = 0.95). The relative amounts and scattering 
characteristics of each population were characterized to quan- 
tify cell morphology. 

Collection of dielectric spectroscopy data 

The culture was monitored by a frequency-scanning per- 
mittivity probe (iBiomass 465 monitor, Fogale Nanotech, 
Nimes, France). The probe was equilibrated in medium 
under relevant culture conditions overnight to ensure no 
probe drift was occurring. A blank medium scan was col- 
lected before inoculation and subtracted from subsequent 
scans, which were collected every 30 s. Subsequent data 
processing was performed in MATLAB (2010a, The Math- 
Works, Natick, MA). Raw scans were filtered for noise using 
a moving average (n = 15 samples). Electrode polarization 
was corrected using the method supplied in the iBiomass 
software. A single-relaxation Cole-Cole model (described in 
the Background section) was fit to the data, which was fil- 
tered and corrected for electrode polarization using the iBio- 
mass software. The dielectric values used for the analysis 
described here are therefore single-term Cole-Cole model 
values supplied by the instrument. 

Prediction of VCV from dielectric increment 

A linear correlation was generated between VCV and per- 
mittivity measured at 1000 kHz. To reduce bias due to shift- 
ing cell states, only samples with viabilities above 95%, as 
measured by trypan blue dye exclusion, were used to build 
the linear correlation. The linear correlation takes the form 
of Eq. 3, where A and k\ are determined as the respective 
slope and intercept of the linear correlation. 

Visualization of multifrequency permittivity spectra 

The dielectric spectra (frequency scans) were analyzed in 
Microsoft Excel 2010 (Microsoft Corp., Redmond, WA) to 
visualize the shape of each frequency scan over time. To bet- 
ter visualize the shape of each scan relative to other scans, 
each frequency scan was normalized from 1 to 0, using 
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where s N (f) is the normalized permittivity as a function of 
frequency and e L is the lowest frequency permittivity mea- 
surement in the scan. 

This method normalizes the dielectric increment, or mag- 
nitude, of each scan, enabling direct comparison of curve 
shapes without the dielectric increment changing due to 
changes in VCV. 

Results and Discussion 

Single-frequency permittivity results 

Single-frequency permittivity values were measured using 
dielectric probes through all phases of a cell culture using 
the GS — /— NS/0 cell line, as described in Materials and 
Methods section. As expected, the permittivity values 
diverged from a linear correlation when plotted against VCV 
during the death phase of the culture as viability declined 
(Figure 2). 

For the growth phase (points above 95% viability), permit- 
tivity measurements exhibited a strong linear correlation 
(R 2 = 0.99) with the measured offline bio volume measure- 
ments, whereas for the death phase (points below 95% via- 
bility) measurements diverged progressively from the linear 
correlation. The divergence manifested as a higher-than- 
expected permittivity value for a given offline VCV mea- 
surement when compared to the growth-phase points. This 
fact supports the initial hypothesis that additional permittiv- 
ity is being measured from a cell population that is excluded 
from the offline viability measurements when viabilities are 
low. An additional published study reports this trend for a 
different cell line (CHO) in a different process from the one 
presented here. 35 

The constants A and k\ in Eq. 3 were determined from the 
linear growth-phase correlation between permittivity and 
VCV for viabilities above 95% to obtain 



VCV=9.17X10" 10 (e)-0.61. 



(6) 



when this equation was used to predict the VCV over the 
course of the entire culture, the death-phase divergence was 
clearly demonstrated, as shown in Figure 2. Any feed-control 
strategy attempting to tailor the feed rate to the amount of 
cell mass present in the reactor would have drastically overfed 
the culture beginning at about 150 h, when the death-phase 
divergence begins to manifest in the permittivity signal. At 
this point in the culture, the product concentration has reached 
only about 60% of the total achieved at the end of the culture 
(data not shown). Any deleterious effects caused by overfeed- 
ing would, therefore, significantly affect the ending titer and, 
possibly, product quality if not corrected. 



Multifrequency permittivity results 

Multifrequency permittivity values were also measured 
during the cell culture run using dielectric probes. To better 
visualize the changing shape of the beta dispersion over the 
course of the culture, frequency scans were normalized from 
1 to 0, as described in the Methods. A sampling of four nor- 
malized frequency scans taken at representative time points 
in the culture highlights the shape change in the scans that 
occurs over the course of the culture (Figure 3). The scans 
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Figure 2. Actual VCV, viability, and predicted VCV values 
using uncorrected permittivities. 
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Figure 3. Normalized frequency scans at representative time 
points during the culture: Day 3 (| |) and Day 6 (O) 
correspond to offline viabilities of >95% (i.e., the 
growth phase). Day 8 (A) and Day 11 (O) represent 
offline viabilities of <95% (i.e., the death phase). 



collected at Days 3 and 6 correspond to the growth phase of 
the culture, when the single-frequency permittivity signal 
correlates linearly with VCV. These scans are so similar in 
shape that they overlap. The scans collected at Days 8 and 
11 correspond to the death phase of the culture, when the 
single-frequency permittivity signal versus VCV diverges 
from the original linear correlation. These scans move pro- 
gressively toward the higher- frequency end of the spectrum 
as the viability of the culture declines. 

The scans collected on Days 3 and 6 correspond to offline 
trypan blue viability measurements of 98% and 95%, respec- 
tively. Scans collected on Days 8 and 11 correspond to off- 
line viability measurements of 70% and 32%, respectively. 

The relative shape of each scan was quantified using a 
novel area ratio (AR) metric from each frequency scan. The 
ratio of area under the portion of each scan above a semiar- 
bitrary frequency, /q, over the area under the entirety of 
each scan is calculated as shown in Eq. 7: 



./// 



AR = 



fa 

~fn 



/h >/q >/l 



(7) 
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where AR is the area ratio for the given scan, f H is the highest 
frequency of the scan, / L is the lowest frequency of the scan, 
and/q is a semiarbitrary frequency chosen between f H and/ L - 

The AR metric was chosen to normalize VCV effects, 
which manifest as a changing dielectric increment, yielding 
equally weighted relative shape information for each scan. 
The AR metric also allows quantification of the relative 
shape of the dielectric spectrum without requiring knowledge 
of the nature of the curve shape a priori. Integrals were cho- 
sen as calculated quantities to decrease the noise associated 
with the method. 

For this work, we chose 808.5 kHz as the /q frequency. 
This was determined to be a good value based on a rough 
sensitivity analysis. Most other chosen frequencies in the 
middle of the scan yielded similar results. The integrals were 
calculated by trapezoidal rule. 

For each offline VCV measurement and corresponding 
single-frequency permittivity measurement, a "required 
correction," or distance from the linear growth-phase model 
defined in Eq. 6, was calculated to determine the amount 
that must be subtracted from each permittivity measurement 
to force all measurements to fall on the linear growth-phase 
line for VCV versus permittivity (Eq. 8): 

^mod — £ — £pred i (8) 

where S mod is the divergence between the linear model 
defined in Eq. 3 and the measured permittivity and e pred is 
the model-predicted permittivity from Eq. 3. This calculation 
requires knowing the offline VCV and using it to build an 
initial calibration. 

Next, <5 mod was predicted based on the AR of each corre- 
sponding frequency scan. The AR from each frequency scan 
(Eq. 7) was linearly correlated to the required signal correc- 
tion (<5 mod; from Eq. 8) to yield 

d mo d(t)=BXAR(t)+k 2 , (9) 

where B and k 2 are the respective slope and intercept con- 
stants generated from the correlation of AR versus S mod . 

Eq. 9 was then fitted to the current data set. The result 
was a linear correlation (R 2 = 0.93), as shown in Figure 4. 
The constants B and k 2 from Eq. 9 were determined from 
the respective slope and intercept of the linear correlation 
between 3 mod and AR (Figure 4) to obtain an expression for 
determining the correction factor from the AR measurement: 

<5mod =78.18XA#-36.04 (10) 

The fact that (5 mod , which quantifies the amount of diver- 
gence observed during the entire culture between the linear 
growth-phase model and the actual measurements, is shown 
to be linearly correlated with the AR metric — a measure of 
the shape of the dielectric spectrum and a direct consequence 
of the dielectric properties of the cell population — supports 
the hypothesis that the divergence during the death phase is 
related to dielectric changes in the cell population. 

An expression for modifying the permittivity (i.e., apply- 
ing the correction factor) for subsequent measurements based 
on the AR was obtained by rearranging Eq. 5 and substitut- 
ing Eq. 9 for S mod to yield Eq. 11: 

e c = 8(t)-(BXAR(t)+k 2 ), (11) 

where s c is the AR-modified permittivity. When the AR- 
modified permittivity is plotted against the offline VCV for 
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Figure 4. AR versus divergence from model for permittivity of 
linear growth phase versus VCV. 
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Figure 5. AR-modified (a) and unmodified (b) permittivity ver- 
sus VCV, showing points above (| |) and below (x) 
95% viability by offline trypan blue dye exclusion. 

the current data set (Figure 5 a), the prediction of all VCV 
measurements, regardless of culture phase, is better explained 
than when using the uncorrected permittivity measurements 
alone (Figure 5b). 

In the analysis presented here, all permittivity measurements, 
regardless of culture viability, are corrected using the method. 
Since the AR values of the frequency scans collected during the 
high-viability growth phase correspond to a very small <5 mod 
value, very little adjustment is applied. The quality of the 
growth-phase prediction for bio volume for points above 95% 
viability (which require no correction), is not significantly 
impacted before (R 2 = 1.00) and after (R 2 = 0.99) analysis. 

To obtain a new expression for prediction of VCV from the 
AR-modified permittivity measurements alone, Eq. 1 1 is substi- 
tuted for s(t) in Eq. 3 to produce the following final expression 

VCV (t) =AX(e(t) -BXAR(t)-k 2 )+h , (12) 

where VCV(f) is the predicted VCV at time t, e(t) is the 
measured permittivity at a single frequency at time t, and 
AR(t) is the measured AR at time t. 
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Figure 6. Predicted and actual VCV values using uncorrected 
and corrected permittivities. 

The measured single-frequency capacitance and AR values 
were then used to produce a corrected VCV using Eq. 12. 
The resulting prediction of VCV vs. time is greatly 
improved, as shown in Figure 6. 

Using the novel methods presented here, VCV can be 
accurately predicted during all phases of the culture, regard- 
less of changing cell health. In addition, this prediction can 
be generated using only data collected from a frequency- 
scanning dielectric probe and an initial calibration curve 
relating A, B, k\, and k 2 to VCV measurements. 

To directly relate changes in the dielectric properties of a 
cell population to the observed changes in the single- 
frequency permittivity measurement of the cell population, a 
single cell-culture data set was used to generate a linear 
model relating the shape of the beta dispersion to the diver- 
gence between actual and predicted VCV measurements. 
This linear model was then applied to the same data set to 
show that a model of this type could account for the diver- 
gence between measured and predicted VCV in the data set. 
Importantly, the application of this method for correcting 
VCV in a bioreactor setting, where prediction of VCV from 
permittivity data is desired based on a previously generated 
calibration, requires that the model constants A, B, k\, and 
k 2 , are constant between runs. Since these constants correlate 
changes in the dielectric properties of the cell population to 
trypan blue dye exclusion, a fundamentally different method 
of determining VCV than dielectric measurements, changes 
in the dielectric properties associated with cell death may be 
expected to change the constants. Therefore, more work 
remains to determine the sensitivity of these constants to dif- 
ferent culture conditions and cell lines. 



Comparison ofAR method to other correction methods 

The analysis presented here depicts one possible method 
for using data obtained from frequency scanning to obtain a 
corrected VCV prediction. Other methods based on multivar- 
iate partial least squares (PLS) modeling have been reported 
by several authors. ' ' Additionally, other methods for 
quantifying the shape of the beta dispersion, such as fitting a 
Cole-Cole model, 16-18 are well established. 

Multivariate PLS methods have the advantage of using the 
latent structure of the data itself to determine the relative 
contributions of input variables (loadings) toward predicting 
the output, rather than defining important variables a priori. 
Because the loadings for multivariate models are determined 
by the latent structure of the data used to build the model, 
the addition of additional datasets typically improves the pre- 
diction of the model. 43 ' 44 This type of approach likely will 
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Figure 7. Actual and predicted VCV values using no correc- 
tion, AR correction, and f c correction methods. 

result in more-accurate predications over other methods 
(e.g., the AR method or Cole-Cole modeling approach) 
when large volumes of data are available, such as late-stage 
process development or manufacturing. 

The AR method is therefore more suitable for use in situa- 
tions where a large calibration data set is unavailable, such 
as in early process development for a cell-culture process. In 
early process development, real-time prediction of VCV is 
often desired for process monitoring and control. 16 

Other techniques for quantifying the shape of the beta dis- 
persion based on fitting Cole-Cole models to dielectric spec- 
tra have been established. To determine whether a similar 
correction using a Cole-Cole model to quantify the shape of 
the beta dispersion could be applied with similar prediction 
capability, the spectra from the current data set were fit to a 
single-term Cole-Cole model. The critical frequency (f c ) 
term of the resulting fit model was used as a surrogate to the 
AR, using the same methods described in the previous sec- 
tion. The results are shown in Figure 7. 

The results of the analysis using f c as the quantifier of 
beta-dispersion shape show that a correction is achieved in 
the advanced death phase. The correction achieved at near- 
peak VCV, as well as in the early growth phase, is consider- 
ably worse than the AR method, however. 

The considerable bias in the early growth phase is likely 
due to poor fitting of the model to the beta dispersion due to 
the low signal obtained from low VCV. This may be 
improved with better model fitting or data-filtering techniques. 

The prediction bias that occurs near the peak VCV is 
likely not attributable to model fitting, because the best 
signal-to-noise ratio occurs at high VCV. The existence of 
this bias suggests that changes in the shape of the beta dis- 
persion, that occur in the early death phase and that are 
quantified by the AR, are not quantified by changes in f c . 
Eventually, at very low viabilities, the f c method does 
quantify shape changes, resulting in a similar correction at 
these points to the AR method. This fact may be related to 
modeling a system that may consist of two discrete relaxa- 
tions as a single relaxation. The prediction may also be 
confounded by the fact that the a term, which also affects 
the shape of the beta dispersion, is also changing (data not 
shown). 



Investigation of cell morphology by flow cytometry 

In addition to implementing the AR method to correct for 
the death-phase divergence, cell samples at different points 
throughout the culture were analyzed via light scattering on 



Biotechnol. Prog., 2014, Vol. 30, No. 2 



485 



a flow cytometer to investigate the biological/physiological 
nature of the link between cell health and the shape of the 
dielectric spectrum. 

Cells were gated into "healthy" and "unhealthy" popula- 
tions according to forward- and side- scattering intensities, as 
described in the Methods section. Histograms of forward- 
scattering intensity and representative visual microscope 
images of cell samples from Day 2 (viability = 97.5%) and 
Day 10 (viability = 32.3%) are shown in Figure 8. Both 
methods show a change in the ratio of healthy to unhealthy 
cells and a distinctly different morphology associated with 
the unhealthy population. This difference in morphology 
between the populations manifests in the roughly Gaussian 
distribution of forward-scattering intensities in the healthy 
population and the clearly bimodal distribution of forward- 
scattering intensities in the unhealthy population. 

Correspondingly, the healthy cells show a roughly uniform 
and circular cell morphology, whereas the size and shape of 
the unhealthy cells is less uniform. The unhealthy cells are 
largely nonspherical and increased cell "granularity" is seen. 
These observations are well-known to be associated with via- 
bility changes in mammalian cells. 41 

Cell morphology changes of this type have previously 
been observed to contribute to shifts in dielectric spectra. It 
has been previously reported 14 ' 45 that cell size has a signifi- 
cant effect on the dielectric increment and the critical fre- 
quency. The observed decrease in cell size would support 
the observed shift of the dielectric spectrum to higher fre- 
quencies. Additionally, the bimodal nature of the forward- 
scattering intensity of the unhealthy cell population suggests 
that the observed shifts in the dielectric spectra may actually 



be a sum of two relatively discrete cell populations. A possi- 
bly fruitful next step could be to explicitly relate the amount 
of healthy and unhealthy cells measured by an offline 
method (such as light scattering) to the dielectric spectrum, 
by splitting the spectrum into two distinct relaxations — one 
relaxation contributed by healthy cells and the other by 
unhealthy cells. This type of analysis could be used to fur- 
ther support the hypothesis that changes in the shape of the 
beta dispersion are caused by permittivity of nonviable cells. 

Changes to overall cell shape that do not involve changes 
in size have also been observed to affect the shape of the 
dielectric spectrum. 46 Elongated, nonspherical cells may also 
explain the observed behavior of the dielectric spectrum. 

In addition to changes in cellular morphology, the link 
between viability and shape of the dielectric spectrum may 
also be explained by intracellular effects. Significant 
changes in the beta dispersion associated with intracellular 
conductivity and membrane capacitance have been 
observed. 25 A shift in the dielectric spectrum associated 
with cell size, lactate consumption, and cell viability has 
also been reported. 14 

The dominant mode of cell death in fed-batch mammalian 
cell cultures is apoptosis. 47 Several investigators have 
observed differences in apoptotic and nonapoptotic popula- 
tions by dielectrophoresis, a technique for measuring the 
frequency-dependent dielectric behavior of single cells. 28 ' 29 
In each study, a marked increase in intracellular conductivity 
accompanied by a decrease in cell size was observed 
between the two populations early in apoptosis. Our group 
has investigated the observability of apoptosis by dielectric 
spectroscopy, which will be published in a separate paper. 
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Figure 8. Forward-scatter intensity histograms gated for healthy (red) and unhealthy (green) cell populations sampled from early in 
the growth phase at Day 2 and late in the death phase at Day 10. 

Representative light microscope images from Days 2 and 10 show distinct morphology changes between the Day 2 (high viability) and Day 10 (low 
viability) samples. 
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Conclusions 

This work demonstrates the relationship between changing 
dielectric properties of a cell population, as quantified by the 
shape of frequency scans, to observed divergences in pre- 
dicted versus actual VCV. A novel AR method for quantify- 
ing the shape of the frequency scans allowed a linear 
relationship to be established between the observed diver- 
gence between predicted and actual VCV and the shape of 
the frequency scans. The work presented here was performed 
on a single data set to demonstrate the relationship. 

The methods presented here may allow use of the full 
dielectric spectrum to predict VCV in bioreactors more accu- 
rately than with single-frequency permittivity methods by 
accounting for changes in cell physiological state that occur 
in the death phase. The ability to measure VCV accurately 
online in bioreactors would enable the implementation of 
process-control schemes (e.g., specific cell-rate feeding) that 
would be more difficult or impossible using other technolo- 
gies. More datasets would be needed, however, to determine 
the robustness of such a method. The relationship between 
changing cell state and measured single-frequency permittivity 
is complex and likely depends on the nature of the cell state 
measured. A comprehensive study across multiple cell lines, 
process conditions, and cell death modes would improve 
understanding of this relationship and the applicability of the 
method to an online VCV monitoring application. 

Additional work also remains to investigate the specific 
physiological changes that are detected by dielectric spectros- 
copy in the death phase. Such fundamental understanding will 
likely promote understanding of other important cellular proc- 
esses that are of interest to the biopharmaceutical industry. 

This initial work also demonstrates the possible utility of 
using commercial frequency- scanning probes for a use 
beyond measuring VCV: the online monitoring of cell state 
in the bioreactor. Where cell physiology is difficult to 
observe — during process-development experiments (e.g., dur- 
ing media and feed development) or in situ in a commercial 
process — dielectric spectroscopy may provide valuable infor- 
mation about cell state (particularly information about cell 
health) that cannot be obtained in near real time. Because of 
the frequency with which the dielectric data may be col- 
lected, this near-real-time feedback of cell-level data also 
may be especially useful for process control. 
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